Phase behavior and interfacial properties 
of nonadditive mixtures of Onsager rods 
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Within a second virial theory, we study bulk phase diagrams as well as the free planar isotropic- 
nematic interface of binary mixtures of nonadditive thin and thick hard rods. For species of the 
same type the excluded volume is determined only by the dimensions of the particles, whereas 
for dissimilar ones it is taken to be larger or smaller than that, giving rise to a nonadditivity 
that can be positive or negative. We argue that such a nonadditivity can result from modelling 
of soft interactions as effective hard-core interactions. The nonadditivity enhances or reduces the 
fractionation at isotropic-nematic (IN) coexistence and may induce or suppress a demixing of the 
high-density nematic phase into two nematic phases of different composition (Ni and N2), depending 
on whether the nonadditivity is positive or negative. The interfacial tension between co-existing 
isotropic and nematic phases show an increase with increasing fractionation at the IN interface, and 
complete wetting of the IN2 interface by the Ni phase upon approach of the triple point coexistence. 
In all explored cases bulk and interfacial properties of the nonadditive mixtures exhibit a striking 
and quite unexpected similarity with the properties of additive mixtures of different diameter ratio. 

PACS numbers: 61.30.Cz, 64.70.Md, 05.70.Np 



I. INTRODUCTION 

In his paper about the isotropic-nematic (IN) transi- 
tion in solutions of monodisperse, rod-like particles that 
interact through a hard, steric repulsion, Onsager briefly 
discussed a possible extension of his results to polydis- 
perse systems^. Since then, a tremendous amount of 
work has been devoted to the study of the influence of 
polydispersity on the phase behavior of such hard-rod 
fluids, both for the case where this polydispersity is of 
the quenched type 2 *^^ and for where it is of the an- 
nealed typeSiLSiS. Focusing on the former, even the sim- 
plest (binary) mixtures consisting of long hard rods that 
differ only in length or diameter exhibit quite nontrivial 
phase diagrams. In addition to the pure isotropic and ne- 
matic phases of various composition and regions of their 
coexistence, the high-density nematic phase can demix 
(and possibly remix) into two nematic phases of different 
composition (denoted Ni and N2). The reason for the 
existence of an IN transition in binary mixtures is the 
same as that in a monodisperse hard-rod fluids, being a 
competition between orientation entropy and entropy of 
packingi*i°u In contrast, the nematic-nematic demixing 
does not involve changes in excluded volume (i.e. packing 
entropy), but rather a competition between orientation 
entropy and entropy of mixingii. Another interesting 
feature is that for sufficiently large size disparity, the two 
distinct nematic phases do not remix even at arbitrary 
high pressure^. 

Unfortunately, it is quite difficult to compare these the- 
oretical findings with results obtained from actual exper- 
iments. Although rod-like particles can be synthesized 
chemically in various waysiS, typically their size distri- 



bution is mono- or bi-disperse only to a first approxima- 
tion. By contrast, suspensions of rod- like viruses such 
as tobacco mosaic virus, M13, pfl and fd are character- 
ized by a high degree of monodispersity, and are there- 
fore attractive model systems, despite the complicating 
factors associated with their fixed physical dimensions, 
their charged nature and the fact that they are not ac- 
tually infinitely rigid but exhibit some degree of bending 
flexibility. 

Recently, however, experimental procedures have been 
developed that allow one to modify the length and the 
diameter of these viruses^, which opens the possibility 
to form binary mixtures of a well-defined bidispersity. In 
particular, one of the methods is based on altering the 
effective diameter of the fd-virus by coating it with the 
polymer polyethylene glycol (PEG). Studies of such bi- 
nary mixtures of thin and thick rods have revealed coexis- 
tence regions of the isotropic and different nematic phases 
(IN2 and INi), as well as a nematic-nematic coexistence 
region (N1N2) and an IN1N2 triple poinliii Although 
some of the gross features of this experimentally deter- 
mined phase diagram are in agreement with theoretical 
predictions based on an extension of Onsager's second- 
virial theory to binary mixtures of hard rods 15-16 , some 
of the experimental and theoretical findings turn out to 
be in sharp contrast with each other. 

According to the theory, mixtures of thin rods (with 
a diameter D\) and thick ones (diameter D2) of equal 
length L should exhibit a spindle-like IN coexistence 
without any nematic-nematic demixing for diameter ra- 
tios d = D2/D1 < 3.8. Experiments, however, point at a 
broad N1N2 coexistence for a diameter ratio as small as 
d ~ 2.0ii. Furthermore, in the interval 3.8 < d < 4.29 



2 



the single nematic phase demixes according to the the- 
ory into two nematic phases N\ and N 2 of different com- 
position, whilst remixing takes place at sufficiently high 
total density (or osmotic pressure), that is, above an up- 
per critical (or consolute) point. Experiments, however, 
reveal a lower consolute point that closes the N\ — N2 
coexistence^, i.e., the N\ — N2 demixing becomes more 
pronounced with increasing osmotic pressure. 

Possible explanations for these differences may well 
be found in the idealisations incurred when modelling 
the virus particles as infinitely elongated, infinitely rigid 
rods that interact with each other only through additive 
hard-core potentials. Indeed, the virus particles are semi- 
flexible and charged, as already alluded to. In addition, 
the grafted polymer coating is soft and hence compress- 
ible, and the length-to-diameter ratio of the rods is at 
best, say, 50. It is important to recall that the second 
virial theory is believed to be exact only in the limit of 
infinite aspect ratios of the rodsi^. The impact of a finite 
length-to-diameter ratio was recently considered within 
an extension of the so-called Parsons-Lee theory to mix- 
tures of hard rodsit This theory does reproduce a lower 
consolute point for mixtures of thin and thick rods albeit 
only if their length (presumed equal) is extremely small. 
A lower consolute point has also been predicted for binary 
mixtures of semi-flexible hard thin rods of unequal thick- 
ness, at least if their persistence lengths and their widths 
do not differ by more than roughly the square root of 
either persistence length over their contour lengthi2*iS. 
However, the predicted isotropic-isotropic demixing is 
not found in the experiments involving the mixtures of 
naked and coated fd virus particleisii The theoryiSiiS 
does anyway not strictly apply to this experimental sys- 
tem, because of the tacit assumption that the length of 
the rods greatly exceeds their persistence length. 

In an attempt to shed light on the issue, we focus on the 
effects that any nonadditivity of the interactions between 
the two kinds of rod might have on their phase behavior. 
Such a nonadditivity emerges naturally if one replaces 
actual soft rod-rod repulsions by effective hard-core re- 
pulsions, characterized by effective hard-core diameters 
that in effect are distances of closest approach. As is 
well known, the screened- Coulomb interactions between 
charged virus particles in an electrolyte solution can be 
reasonably approximated by an effective steric interac- 
tion with a hard-core diameter that is the sum of the 
bare, "physical" diameter of the rod and an electrostatic 
contribution proportional to the Debye screening length 
of the suspending medium^. For the interaction between 
a pair of polymer-coated virus particles one would have 
an effective diameter of the order of the radius of gyra- 
tion of the tethered chainsSI, at least if the Debye length 
is much smaller than that. 

It is not at all obvious that the effective interaction 
length between a bare and a polymer-coated charged rod 
should be the linear average of the interaction lengths of 
the two separate species, in other words, one would from 
the outset expect the interaction within such an effec- 



tive description to be non-additive rather than additive. 
Indeed, as we shall see below in section even highly 
simplified model potentials produce non-additive effec- 
tive hard-core interactions in mixtures of rods. The level 
of nonadditivity may be expressed in a parameter a de- 
fined such that the effective hard-core diameter of an un- 
like pair of rods can written as \{D\ +Z?2)(1 + &), where 
D a is the effective hard-core diameter of the interaction 
between two like rods of species a — 1,2. For an addi- 
tive mixute, a = 0. In this paper, we make plausible by 
explicitly considering the steric interactions between the 
various types of rod that, even within a simplified model, 
a may attain values that can be positive or negative up 
to, say, ten per cent. Additional sources of nonadditiv- 
ity may be found, say, in electric polarisation effects of 
the charges on the polymer coating, but these will not be 
considered here. 

Although the microscopic origin of nonadditivity is ul- 
timately based on the less (a > 0) or more (a < 0) 
efficient packing of the mixture compared to the pure 
species, we do not attempt to calculate a from a realis- 
tic microscopic theory. Having ascertained that a need 
indeed not be zero, we treat it as a phenomenological pa- 
rameter in a generalised Onsager theory, and investigate 
its consequences for the phase behavior of the mixture, 
and for the interfacial properties of co-existing isotropic 
and nematic phases. As we shall see, both the predicted 
phase diagrams and interfacial properties of the isotropic- 
nematic interface are very sensitive to values of \a\ as 
small as a few per cent. Of course, this does not im- 
ply that all effects of electrostatic interactions, flexibil- 
ity, etc., are accurately or even properly modelled. In 
fact, we find that non-additivity cannot explain the exis- 
tence of the lower consolute point found by Fraden and 
co-workersia. 

The remainder of this paper is organized as follows. In 
Sec. [n]we introduce a simple model for polymer-coated 
rods, and provide an estimate for typical values of the 
nonadditivity parameter a. In Sec. IIIII we introduce 
the Onsager-type free energy functional, and derive from 
that the basic Euler-Lagrange equations describing the 
orientational and density distribution of the rods under 
conditions of thermodynamic equilibrium. In Sec. II VI we 
solve these equations for bulk geometries, and analyze the 
structure of a few typical bulk phase diagrams. In Sec. |V] 
we briefly describe a method to solve the Euler-Lagrange 
equation for interface geometries of binary mixtures, and 
study IN±, N1N2 and IN2 interfaces, the latter in par- 
ticular in the vicinity of the bulk IN1N2 triple point. A 
summary and discussion of the results are presented in 
Sec. ED 



II. NONADDITIVITY OF INTERACTIONS 

The physical origin of nonadditivity can be illustrated 
on the basis of a simple model for a mixture of bare and 
PEG-coated fd viruses^. The bare rods are modelled 
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as rigid hard rods of length L and diameter Ai (L ^> 
Ai), and hence the interaction potential between two 
bare rods, 0n(r), is given by P4>n(r) — oo, for r < Ai 
and r > Ai, respectively. Here, r denotes the shortest 
distance between the main axes of the two rods. 

The PEG-coated rods are identical to the bare ones, 
except that they bear an additional soft layer extending 
to a distance A2/2 from the axis of the rod, i.e., to a 
distance (A2 — Ai)/2 from their hard-core surface. We 
do not specify the relation between the dimensions of the 
tethered PEGs and A2 in any detail, but one expects that 
(A2 — Ai)/2 is of the order of the radius of gyration of 
the grafted PEG (so we only consider A 2 > Ai). We ex- 
pect that the soft, repulsive interaction that occurs when 
the polymer coating of two rods overlap should be quite 
similar to that of overlapping star polymers 2 ^. In order 
to keep the model as simple as possible, we represent the 
interaction of mean force resulting from the presence of a 
polymer coating by a square-shoulder potential that is a 
function of r alone, and ignore any angle dependence that 
might arise in reality. This angular dependence should 
be significant only for configurations of rods inclined at 
small angles, which bear only a tiny statistical weight 
in the limit of large aspect ratios. Note that although 
our representation of the soft potential is isotropic, the 
virials based on it are anisotropic because the interaction 
volumes are a function of the relative orientations of the 
rods. 

The interaction potential between a bare and a coated 
rod, </>i2(r), should obviously be identical to the naked- 
rod potential </>n(r) if r < Aj and r > (A 2 + Ai)/2. 
Within our description, 4>i2 (r) takes on a value different 
from that, £1 > 0, if Ai < r < (Ai + A 2 )/2, i.e., when 
the hard-core of the bare rod perturbes the soft outer 
layer of the coated rod. It is to be seen as an average 
of the actual interaction potential over its range. Our 
effective interaction potential between two coated rods, 
4>22 (r) , is more complicated and consists of two shoulders 
in between the range of the hard-core repulsion (r < Ai) 
and the noninteracting long-distance regime (r > A 2 ). 
The first shoulder, for Ai < r < (Ai + A 2 )/2, is such 
that (j>22( r ) = 2ei, and represents the overlap of the hard 
core of the first rod with the polymer layer of the second 
one, and vice versa by symmetry. The second shoulder, 
for (Ai + A 2 )/2 < r < A2, represents overlap of the two 
polymer layers, and is such that 4>22{r) = £2 > 0. 

The nature of the polymer chains is such that we ex- 
pect their entropy to be reduced more by a penetrating 
rigid rod than by another polymer. Indeed, the cross 
virial of a rod and a flexible chain is much larger than 
the geometric average of the rod-rod and the chain-chain 
virials 2 —. For this reason we only consider cases where 
2ei > e 2 . The pair potentials 4> a a'{ r ) between rods of 
species a and a' are illustrated graphically in Fig. 

It is a straightforward exercise to calculate the second 
virial coefficients B aa i averaged over all angles, from the 
pair interactions 4> aa i (r) given abov o 1 ^ . In the Onsager 
limit L > A2 > Ai, where terms of order LA 2 can be 



ignored, one finds 
B u = (7r/4)£ a Ai, 



B u = (vr/4)^ Ax 



B22 = Ax 



A2-A1 



A 2 -Ai 



1-e 



-/3ei 



(1) 

(2) 



'2 - e- 2l3ei - e~^ 



(3) 



These expressions can be used to map the model mixture 
of bare and PEG-coated rods onto a mixture of hard rods 
with effective hard-core diameters D\ and D 2 . We choose 
D\ and D2 to be such that the like-like second virial 
coefficients of the effective hard-core system are identical 
to Bu and B22 given in Eq. and respectively, 
i.e., we impose that B aa — (Ti/A)L 2 D r ,. This yields 

D x = Ai, 

D 2 = Ai + A " 2 Al (2 - exp(-2/? £l ) - exp(-/3e 2 )). 

(4) 

One may verify that D2 = A2 in the limit that /3ej — » 00, 
as expected. We now also impose that the cross virial 
coefficient of the effective hard-core system equals B\2 
given in Eq. @. For arbitrary t\ and £2 this re- 
quires a nonadditivity parameter a such that B12 = 
(tt/8)L 2 (D 1 +D 2 )(1 + a), which yields 



f 



2(rf-l)(l-exp(-/?ei)) 
2 - exp(-2/? ei ) - exp(-/3e 2 ) 



1, (5) 



with d = D2/D1 the effective diameter ratio of the rods. 
In Fig. [21 we show the contour plot of the nonadditiv- 
ity parameter a as a function of the energy scales £1 
and e 2 for the effective diameter ratio d = 3.5, which 
is the value that we will use in our calculations below; 
other values for d produce similar contour plots. The 
grey area in Fig. [3 is the regime deemed unphysical, 
with 2ej < £2- As one can see, for physically reason- 
able values of e\ and €2 of the order of ksT both pos- 
itive and negative values for a are possible, even when 
£1 > £2- The crossover from positive to negative nonad- 
ditivity takes place, independently from the value of d, 
when exp(— (3e 2 ) — 2exp(— f3e{) — exp(— 2/3ei), i.e., when 
/3£ 2 - (/3£i) 2 if (3ei < 1 and f3e 2 ~ /3ei - In 2 if pei > 1. 
Presuming that both (3e\ and (3e 2 are indeed of the order 
unity, we expect \a\ to be in the range 1CP 2 — 1CP 1 . Such 
small deviations from additivity are sufficient to qualita- 
tively alter the phase behavior of the rods, as we shall 
see next. In our study, we from now on treat a, D\ and 
Z?2 as independent parameters. We investigate both the 
bulk and the interfacial behavior of the effectively purely 
hard-core system, in which the soft interactions are in- 
corporated through the degree of non-additivity a. 
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III. DENSITY FUNCTIONAL AND METHOD 

Consider a fluid of hard cylinders of two different 
species a = 1,2 of diameter D a and equal length L 
{D a /L — > 0) in a macroscopic volume V at temperature 
T and chemical potentials p a . Let r denote the center- 
of-mass coordinate of a rod and u) the orientation of the 
long axis. The interactions between the cc'-pair of rods 
with coordinates q = {r, ui} and q 1 = {r',cD'} are char- 
acterized by a hard-core potential, which is the simple 
contact potential for rods of the same species (er = a'), 
whereas for unlike rods (a ^ a') it corresponds to in- 
teractions between hard rods of diameter (1 + a)Di and 
(l + a)D 2 . 

Within the second virial approximation and in the ab- 
sence of external potentials, the grand potential func- 
tional f^l/v}] of the one-particle distribution functions 
p a (r,Lu) can be writte n 1 ; 10 ! 16 as 

f3Q[{pA] = f dq P M){HpM)L 2 D a ]-l-P^ 

J d qdq' f <J a'{q\q')Pa{q)Pa'{q'), ( 6 ) 

with (3 = (fcsT) -1 the inverse temperature, and 
fuu'{q,q') the Mayer function, which equals —1 if the 
rods overlap and vanishes otherwise. Since we consider 
the limit D a / L — > for any a, the relative shape dispar- 
ity of rods is characterized by the ratio d = D 2 jD\ of the 
diameters and the value of the nonadditivity a. 

The minimum conditions 5VL[{p a }] / 5 p a {q) = on the 
functional lead to the set of nonlinear integral equations 

ln[p„(q)LlD„]-J2 f dq' UAvrfWW) = P»o (7) 

cr' 

to be solved for the equilibrium distributions p a {q). 
These equations are identical to the Euler-Lagrange 
equations for additive rods mixtures, and we can directly 
apply the method developed earlieri^i 6 .. The structure 
of the bulk phase diagram depends now on the value of 
the nonadditivity parameter a, and has to be determined 
first. 

Since the bulk distribution functions of the isotropic 
and nematic phase are translationally invariant, i.e., 
p a (r,Lu) = we can reduce Eq. (7J to 

hx\p a (Cj)I? a D a ] J du'E aa , (&,&') p<t>(&') = Ppcr, 

a' 

with E aa i the excluded volume of a pair of cylinders of 
species a and a 1 given by 

E aa ,{u,Cb') = - J dr'f aa/ (r,u);T',u)') 

= L 2 (Ar + ArO(l + a(l-<W))|sin<^| 

(9) 



in terms of the angle if between Qj and uj' ', i.e., ip = 
arccos(ti) • ui'). Note that additional 0(LD 2 ) terms are 
being ignored in Eq. Q, in line with the needle limit 
{D a /L ~ > 0) of interest here. Given the linear depen- 
dence of the excluded volume on D a , one can see that 

E 12 {lu,lu') = ^(E n (u>,oj') + E 22 (u,u')) (1+a). (10) 

In some sense, a plays a similar role in the present context 
as the so-called \ -parameter in the Flory theory of poly- 
mer solutions on a lattice, where demixing is driven by 
direct unfavorable nearest neighbor interaction between 
unlike species as compared to that between like species. 

Details of the numerical schemes to solve Eq. (jHJl have 
been discussed elsewher o 15 ! 16 . Here we use a nonequidis- 
tant 6>-grid of Ng — 30 points B. L e [0, tt/2], where 
1 < i < Ng, in order to find the bulk distributions p a (9i). 
Coexistence of different phases {/, N 1} N 2 } is determined 
by imposing conditions of mechanical and chemical equi- 
librium. 



IV. BULK PHASE DIAGRAMS 

In Fig. |3 we show both pressure-composition (a) and 
density-density (b) representations of bulk phase dia- 
grams of thin-thick binary mixtures (L a = L, D 2 > D\) 
for the diameter ratio d = 3.5 at various values of the 
nonadditivity parameter a. In Fig. EUa) the composi- 
tion variable x = n 2 /(n\ + n 2 ) denotes the mole fraction 
of thick rods, n a — J du>p a {&) is the number density 
of species a, and p* — (tt / '4)(3pL 2 D\ is a dimension- 
less bulk pressure. Note that the IN coexistence pres- 
sure pthin and pthick of the pure thin {x = 0) and pure 
thick (x = 1) system are given by (ir/4)(3p t hinL 2 Di — 
(Tr/4)f3pthickL 2 D 2 = 14.045, i.e., pthick = Pthin/ d, and 
that the tie-lines connecting coexisting phases are hori- 
zontal in the p — x representation of Fig. I^a). This rep- 
resentation is convenient for theoretical analysis, whereas 
the densities (volume fractions) of thin and thick rods are 
experimental control parameters^. For this reason the 
same phase diagrams of thin-thick binary mixtures are 
shown in Fig. CUb) in the density-density representation, 
with nj = niL 2 Di(ii/A) and n* 2 = niL 2 D 2 (n/A) being 
the dimensionless bulk number densities of thin and thick 
rods, respectively. In this representation the tie-lines, in- 
dicated by the dotted lines, are no longer horizontal. 

The structures of the bulk phase diagrams for various 
a show a striking similarity with the bulk phase diagrams 
of additive binary mixtures of thin and thick hard rods 16 . 
At low pressures (or low densities) the phase diagrams 
show an isotropic (J) phase, and at higher pressures (or 
densities) one (a < 0.07) or two (a > 0.07) nematic 
phases (TVi and N 2 ). For a < 0.07 the phase diagram 
is spindle-like, and the only feature is a strong fraction- 
ation at coexistence, such that the nematic phase is rel- 
atively rich in thick rods and the isotropic phase in thin 
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ones. Although the nonadditivity modifies the fraction- 
ation gap, the reason behind it remains the same: the 
relatively large excluded volume in interactions of the 
thick rods makes them more susceptible to orientational 
orderingi 5 *^ 5 ^^. As a general tendency, the fractionation 
at isotropic-nematic coexistence becomes stronger for in- 
creasing values of a. 

For a > 0.06 the bulk phase diagram develops nematic- 
nematic (N1N2) coexistence in a pressure regime p > pt, 
with pt the triple-point pressure. Using the simple Gaus- 
sian ansatz for one-particle distribution functions, one 
can demonstrate that the packing entropy does not play 
a role in nematic demixing in our system, similar to the 
case of additive mixturesii. Although it is known that 
the functional form of p a (&) is not Gaussian even at high 
densities, an analysis of the exact high-density distribu- 
tion functions confirmed such a mechanism of nematic 
demixing 15 . On this basis we assume it to be valid at 
arbitrary high pressure in our system, and expect the 
structure of the bulk phase diagrams to be similar to 
those of additive mixtures. In particular, for a — 0.07 
nematic remixing is observed at a sufficiently high pres- 
sure, as illustrated in Fig. The consolute point, at 
which the density and composition difference between the 
coexisting nematic phases vanishes, is indicated by (*). 
For a = 0.1, limitations of the numerical scheme^! do 
not allow us to determine whether or not remixing takes 
place at high enough pressures. We note that in the limit 
of very high pressures, where the rods increasingly align 
themselves, both end corrections and higher order virials 
need to be taken into account for an accurate description 
of the phase behavior. On the other hand, in analogy 
with additive mixtures, one expects that critical values 
of a and d, beyond which the nematic demixing does not 
take place at arbitrary high densities, existi 5 -. 

In order to characterize the amount of nonadditivity in 
the excluded volume interactions which leads to signifi- 
cant structural modification of the phase diagram (i.e. 
nematic demixing), we explore various thin-thick mix- 
tures of different values of a and d, and determine the 
value of a* for which the pressure of the nematic- nematic 
consolute point and the triple point pressure coincide. 
Results of our studies are presented in Fig. 0] For a < a* 
(at fixed d) the N1N2 phase separation is not detected, 
and for a > a* there is N1N2 coexistence in the phase 
diagram. It is evident that in the interval d € [3.5,4.2] 
even a small nonadditivity \a\ < 5 — 7% may induce or 
suppress the N1N2 demixing transition. 

One might surmise that the linearity of the function 
a*(d) within the explored range of a reflects the linear- 
ity of the excluded volume Eyi (u), Cj') in terms of d and a, 
because it drives the nematic-nematic phase separation. 
The mapping of the non-additive to the additive case is 
not trivial, however, since the density distributions p a (<2)) 
depend on a implicitcly. Nonetheless, direct comparison 
of the bulk phase diagrams of the nonadditive mixture 
with d = 3.5 and a = 0.07 and the additive mixture 
with, for instance, d = 4.0i& shows close values of the 



fractionation gap at the N1N2 coexistence. Further evi- 
dence for similarity of these systems in the high density 
regime will be demonstrated in our analysis of their in- 
tcrfacial properties presented next. 

V. INTERFACES 

Free planar interfaces between various coexisting bulk 
phases can be studied similar to the interfaces of addi- 
tive mixturesiSiSSiS. We focus on the nonadditive thin- 
thick mixture characterized by d — 3.5 and a = 0.07. 
The nematic director n of the asymptotic nematic bulk 
phase(s) can, in general, have a nontrivial tilt angle 
8 t = arccos(n • z) with respect to the interface normal 
z. In the present calculations we restrict attention to 
9 t = 7r/2, i.e., n _L z. As we have checked, this geometry 
is thermodynamically favorable because of its minimal 
surface tension. 

Similar to the studies of additive mixtures, we use the 
planar symmetry of the interfaces and assume the distri- 
bution functions to be uniaxially symmetric with respect 
to the director, i.e. /? CT (r,(£>) = p a {z,9), which reduce 
Eqs. to 

f3[x a = \n[p a (z,8)LlD a }+Y^ J dz'd6'sm8' 

a' 

x/Cw(z-z',M'W(^'), (11) 

with K a<T i{z - z',9,9') = / difdtp'dx'dy' f a<T >(q,q'). 
We solve Eq. in order to determine uniaxially 

symmetric nonuniform distributions p a (z,9i) using an 
equidistant z-grid of N z = 200 points in the inter- 
val z G [— 5L, 5L], and corresponding bulk distributions 
p a {9i) as boundary conditions. Further details of the nu- 
merical calculations were discussed elsewhere^. 

The IN\ and N1N2 interfaces are found to be smooth 
and monotonic, in the sense that the profiles of the ne- 
matic uniaxial order parameters S a (z) and the densities 
n a (z) change monotonically from the bulk values in the 
/ (Ni) phase to those in the N% (N2) phase. The cor- 
relation length of the bulk iVi phase at the triple- 
phase coexistence (as well as £/ and £/v 2 for the / phase 
and the N% phase, respectively) can be extracted from 
the asymptotic decay of the one-particle distributions 
p a (z,6) to their bulk values p^ioj), since the deviation 
5p a {z,Cj) = p a (z,Lu) — p^ x {uj) is of the formic 

5p a (z,tb) = A T (o))exp(-z/£iv 1 ), z ~> 00 ■ ( i2 ) 

Interestingly, we find that £,n x /L = 0.49 ± 0.02 is the 
same as for the additive mixture with d = 4.0 16 , which 
has a virtually identical phase diagram. 

The properties of the IN 2 interfaces depend strongly 
on the pressure difference with the triple point (IN1N2 
phase coexistence). As it is demonstrated in Fig. 
the surface tension of the IN2 interface shows a non- 
monotonic dependence on the bulk pressure p, and 
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strongly correlated with the fractionation at the IN2 co- 
existence. Upon increasing the nonadditivity, the surface 
tension 7jtv 2 (p) grows, again indicating that a plays a 
role similar to the diameter ratio d. For comparison we 
have included 7/jy (p) for an additive thin-thick mixture 
with d = 4.0, which is again quite close to the results for 
the nonadditive mixtures with d = 3.5 and a — 0.07 and 
0.1. 

The microscopic thickness t of the interface is defined 
as t = \z + — Z—\ where z± are solutions of n\"(z) = 0, 
and a prime denotes differentiation with respect to z. 
As this equation has a set of solutions in every interfa- 
cial region, we choose for z± the outermost ones, i.e., 
the ones nearest to the bulk phases. The density of thin 
rods is a convenient representation of structural changes 
within the interface, since they have a smaller excluded 
volume and a non-vanishing concentration in both coex- 
isting phases. This criterion provides a single measure 
for the thickness of both monotonic and non-monotonic 
profiles, with and without a thick film in between the 
asymptotic bulk phases at z — > ±00. The interfacial 
width for the one-component IN interface is, with the 
present definition, given by t/L = 0.697. 

The thickness of the IN% interface was found to diverge 
upon approach of the triple-point pressure pt ■ This can 
be seen in Fig. H3 where t/L is plotted as a function of 
the dimensionless undersaturation e = 1 — p/pt, which 
is a convenient measure of the pressure difference with 
the triple point. The nature of the film can be analyzed 
from the density profiles 711(2) of the IN2 interface (or 
equivalently S a (z), or 712(2))- ln Fig. 0the profiles of 
711(2) are shown at several values of the undersaturation 
e. The asymptotic densities at z — > ±00 in Fig. [7| arc 
those of the coexisting / and N2 bulk phases (at the cor- 
responding e). For e — > the value of 711(2) in the film 
approaches the density of thin rods of the bulk triple 
point N\ phase, indicated by the dashed line in Fig. 
However, the undersaturation e = 10~ 4 is yet too large to 
be in the asymptotic thick-film regime. The same iden- 
tification can be made for 712(2) and S a (z), and on this 
basis we conclude that the wetting phenomenon under 
consideration is complete triple point wetting of the free 
IN2 interface by the N\ phase. The similarity with com- 
plete wetting of the IN2 interface by the N± phase in 
additive thin-thick mixture with d = 4.0 is again rather 
striking, as is clear from Fig. El where the thickness of 
the IN2 interface for additive rods is indicated by (o). 

Since one expects, for the short-range interactions of 
interest here, that the thickness of the wetting Ni film in 
the IN2 interface diverges as t ~ —^N! hie for e — > O^S, 
the value of the correlation length of the bulk N± phase 
can be extracted, £/Vi = 0.49 ± 0.02, which is consis- 
tent with the value determined earlier from the decay of 
p a (z, 9) into the bulk N± phase. 

The analysis of the structural properties of the IN2 
interface can be complemented by studies of the ratio of 



surface tensions 



lir %lp 4 (7/JVi +7iViiV 2 )' 

as presented in the Inset in Fig. It is clear that upon 
approach of triple-point coexistence lim c ^o R( e ) = lj 
which implies a vanishing contact angle. This provides 
the thermodynamic proof of complete triple-point wet- 
ting. The e-dependence of R again reveals that e = 10~ 4 
is too large to be in the asymptotic thick-film regime. 

VI. SUMMARY AND DISCUSSION 

In this paper we have explored the bulk phase dia- 
grams and the interfacial properties of the nonadditive 
mixtures of thin and thick hard rods. The nonadditiv- 
ity was introduced in an attempt to effectively capture 
some of the effects of soft interactions between them, hav- 
ing in mind mixtures of bare and PEG-coated fd virus 
particles in aqueous suspension 1 ^. We showed that the 
effective hard-core diameter of the unlike interactions, 
\{D\ + D 2 )(l + a) with D\ and D 2 the effective diame- 
ter of the like interactions, can easily be smaller or larger 
than the additive case, |(-Dj + D 2 ), by more than a few 
per cent. 

As is illustrated in Fig. |3 a small amount of non- 
additivity a > can stabilize the high-density nematic- 
nematic phase coexistence, even if it is only metastable 
for an additive mixture with the same diameter ratio. 
However, the experimentally observed lower critical point 
of the nematic-nematic demixing transition 1 ^ could not 
be reproduced by incorporating nonadditivity into the 
theory. We suggest, therefore, that further theoretical 
studies of this system should consider in more detail the 
impact in particular of a finite bending flexibility, be- 
yond the ground-state approximation 1 ^. Another issue 
that needs to be resolved is the (elastic) response of the 
polymer coat to volume exclusion between the rods, an 
aspect completely ignored in our analysis. 

We present results of bulk and interfacial calculations 
for the specific diameter ratio d — 3.5, motivated by the 
experimental parameters 1 ^. Other values of d lead to 
similar conclusions. We find the bulk phase diagrams 
of nonadditive binary mixtures to show a large similarity 
with those of the additive mixtures of larger diameter ra- 
tio. This is most likely related to the linear dependence of 
the rod-rod excluded volume on both the diameter ratio 
d and the nonadditivity a although it is not clear whether 
there is an exact mapping linking non-additive and ad- 
ditive hard-rod mixtures. We also found that many if 
not all of the interfacial phenomena that we studied are 
similar to those of additive mixtures with a larger di- 
ameter ratio. Similar to the interfaces between different 
bulk phases in additive mixtures, the IN\ and N1N2 in- 
terfaces are smooth and monotonic, whereas the IN2 in- 
terface exhibits complete wetting by the Ni phase upon 
approach of the triple phase coexistence. The complete 
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triple-point wetting scenario was confirmed by (i) the log- 
arithmic divergence of the thickness of the N± film with 
vanishing undersaturation, and (ii) the surface tension 
ratio lim e ^o R= 1- Such a similarity between properties 
of additive and nonadditive mixtures may represent a 
significant difficulty to distinguish these in experiments. 
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FIGURE CAPTIONS 

1. Interaction potentials <j>u(r), (f>i 2 (r), 4>22{r) between bare-bare (solid), bare-coated (dashed), coated-coated 
(dotted line) rods, respectively, as a function of the (shortest) distance r between the axes of the two rods. 

2. Contour plot of the nonadditivity parameter a as a function of the square-shoulder values t\ and e 2 for the 
effective diameter ratio D 2 /Di = d = 3.5. The grey area denotes the nonphysical region e 2 > 1t\ (see text). 

3. (a) Bulk phase diagrams of binary thin-thick mixtures (diameter ratio d = 3.5) for different nonadditivity param- 
eter a in the p — x representation, with p* — (ir/4)(3pL 2 Di the dimensionless pressure, and x the mole fraction 
of the thicker rods. We distinguish the fully symmetric isotropic phase (J) and orientationally ordered nematic 
phases (N\ and N 2 ). For the nonadditivity parameter a — 0.07 the IN1N2 triple phase coexistence is marked 
by (A), and the NiN 2 critical point by (*). (b) The same phase diagrams in density-density representation, 
where n\ — n\L 2 Di(n j '4) and n 2 = n 2 L 2 D 2 (ir / 4) are the dimensionless bulk number densities of thin and thick 
rods, respectively. The tie-lines connect coexisting state points. 

4. Nonadditivity parameter a* at which the consolute point and triple point coincide for various values of the 
diameter ratio d. For mixtures, characterized by a(d) < a*(d) the N\N 2 demixing is not detected. 

5. Dimensionless surface tension ^} N2 = Pjin 2 LD\ of IN 2 interfaces as a function of dimensionless pressure 
p* = f3pL 2 Di(ir/4) for thin-thick mixtures with diameter ratio d = 3.5 and different values of nonadditivity 
a = 0.0 (o), 0.07 (o) and 0.1(D). The dashed line indicates j} N2 (p*) for the additive thin-thick mixture of 
diameter ratio d = 4.0. 

6. Thickness t/L as a function of the undersaturation e = 1 — p/pt from the triple point pressure pt for thin-thick 
mixtures with diameter ratio d = 3.5 and nonadditivity a — 0.07 (*). For comparison we show thickness of the 
/./V2 interface of the additive mixtures for d = 4.0 (o). The inset shows the surface tension ratio R [see Eq. Q13|l] 
as a function of the triple point undersaturation e. 

7. Density profiles of the thin rods n*(z) in the IN 2 interface for diameter ratio d = 3.5 and nonadditivity a = 0.07 
at triple point undersaturations e = 1 — p/pt = 10~ 2 , 1Q -2 ' 5 , 10 -3 , 10 -3 5 , 10~ 4 (from bottom to the top curve). 
The bulk I/N 2 phase is at z — > — 00/00. The dashed line n\ = 3.977 represents the bulk density of thin rods in 
the triple point N\ phase. These profiles indicate the formation of a wetting N\ film in the IN 2 interface. 
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